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ABSTRACT 

One of the factors influencing the spectral evolution of a planetary nebula is the 
fate of the dust grains that are emitting the infrared continuum. Several processes 
have been proposed that either destroy the grains or remove them from the ionized 
region. To test whether these processes are effective, we study new infrared spectra of 
the evolved nebula NGC 6445. These data show that the thermal emission from the 
grains is very cool and has a low flux compared to H/3. A model of the ionized region is 
constructed, using the photo-ionization code CLOUDY 90.05. Based on this model, we 
show from depletions in the gas phase elements that little grain destruction can have 
occurred in the ionized region of NGC 6445. We also argue that dust-gas separation in 
the nebula is not plausible. The most likely conclusion is that grains are residing inside 
the ionized region of NGC 6445 and that the low temperature and flux of the grain 
emission are caused by the low luminosity of the central star and the low optical depth 
of the grains. This implies that the bulk of the silicon-bearing grains in this nebula 
were able to survive exposure to hard UV photons for at least several thousands of 
years, contradicting previously published results. 

A comparison between optical and infrared diagnostic line ratios gives a marginal 
indication for the presence of a t^-effect in the nebula. However, the evidence is 
not convincing and the differences could also be explained by uncertainties in the 
absolute flux calibration of the spectra, the aperture corrections that have been applied 
or the collisional cross sections. The photo-ionization model allowed an accurate 
determination of the central star temperature based on model atmospheres. The 
resulting value of 184 kK is in good agreement with the average of all published 
Zanstra temperatures based on a blackbody approximation. 

The off-source spectrum taken with LWS clearly shows the presence of a warm 
cirrus component with a temperature of 24 K as well as a very cold component with 
a temperature of 7 K. Since our observation encompasses only a small region of the 
sky, it is not clear how extended the 7 K component is and whether it contributed 
significantly to the FIRAS spectrum taken by COBE. Because our line of sight is in 
the galactic plane, the very cold component could be a starless core. 



Subject headings: stars: AGB and post-AGB — circumstellar matter — stars: 
fundamental parameters — planetary nebulae: individual (NGC 6445) — infrared: 
ISM: continuum 



1. Introduction 

When low- and intermediate mass stars are approaciiing the end of their evolution, they go 
through a period of heavy mass loss known as the Asymptotic Giant Branch (AGB) stage. In 
this phase the star is a red giant which consists of a degenerate C/0 core surrounded by a very 
tenuous envelope. During the AGB stage the star loses a substantial fraction of its initial mass, 
eventually leaving only the degenerate core of typically 0.6 Mq to 0.8 Mq covered with a thin 
layer of hydrogen. This core will eventually evolve into a white dwarf. It is generally assumed 
that dust formation is very efficient in the AGB outflow so that large amounts of dust grains are 
formed during that phase. Once the envelope mass of the central star drops below a critical value, 
the mass-loss rate will drop several orders of magnitude and the star will enter the post-AGB 
phase of its evolution. During this phase the central star will heat up and the AGB ejecta will 
expand, forming a detached shell. Once the central star reaches a temperature of 20 kK to 30 kK, 
it will start to ionize the AGB shell and a Planetary Nebula (PN) will form. After this the central 
star will continue heating up at roughly constant luminosity until finally the nuclear reactions 
cease and the central star will decrease both in temperature and luminosity. All this time the 
detached AGB shell will continue to expand and undergo hydrodynamic interactions with the fast 
wind emanating from the central star. When the central star enters the cooling track, the number 
of ionizing photons emitted by the star will decline sharply. If this decline in luminosity is fast 
enough (i.e., if the central star is massive enough), this will cause the ionized region to shrink 
and the outermost parts of the nebula will start recombining. Eventually the cooling rate of the 
central star will slow down and the nebula will dissolve in the Inter-Stellar Medium (ISM), leaving 
a bare white dwarf. Both the evolution of the central star and the expansion of the AGB ejecta 
will have a profound effect on the infrared emission from the dust grains, causing a substantial 
evolution of both the total infrared flux and the temperature of the dust grains. 

During the post-AGB phase, the combined effects of increasing stellar temperature (which will 
increase the grain temperature) and nebular expansion (which will decrease the grain temperature) 
cannot be predicted easily and detailed modeling is necessary. A comprehensive discussion of the 
evolution of the infrared spectrum during this phase can be found in van Hoof, Oudmaijer & 
Waters (1997a). However, after the central star has reached a temperature of 40 kK to 50 kK, 
the effect of the rising stellar temperature on the grains will be reversed. This is caused by the 
fact that the peak of the stellar spectrum shifts to ever shorter wavelengths, beyond the peak of 
the grain absorption cross section, thereby making grain heating less efficient. Therefore, once the 
nebula is in the PN phase, the evolution of the infrared spectrum is clear. Both the ratio of the 
total infrared flux over the H/3 flux and the average dust temperature will decrease with increasing 
(Stromgren) radius of the nebula. However, a considerable amount of scatter around this average 
evolution can be expected in a statistical sample, due to the particulars of the evolution of each 
individual nebula. 

One of the factors influencing the evolution of the infrared spectrum is the dust grains 
themselves. In the previous discussion we implicitly assumed that nothing much happened to 



the grains; however, this is not at all clear. It has been proposed that the dust grains could be 
destroyed as the PN evolves (e.g., Natta &: Panagia 1981; Pottasch 1987), either due to spallation 
by hard UV photons or shocks. It is known that UV photons are abundantly present in PNe, as is 
the case for shocks which are formed due to the interaction of the fast wind with the slow AGB 
ejecta. Hence, if these processes are efficient, the ionized regions of evolved PNe should be free of 
dust grains. Another possibility is that dust grains are separated from the gas due to radiation 
pressure. This could lead to a general flow of the grains with respect to the gas or, if instabilities 
occur, to the formation of small, high-density globules which would be neutral or partially ionized. 
Such globules have been observed, suggesting that at least in some nebulae this mechanism is 
effective. 

To investigate these points in more detail, we shall analyze the infrared spectrum of an 
evolved PN and search for indications of any of the processes we mentioned above. This nebula 
is NGC 6445 and it is classified as a bipolar type I PN (Corradi & Schwarz 1995; Perinotto 1991 
respectively), a class assumed to originate from massive progenitors. Based on its low density, the 
nebula appears well evolved. Both the ionized and neutral material (CO) show a high expansion 
velocity of 38 kms~^ (Weinberger 1989) and 33 kms~^ (Huggins & Healy 1989) respectively. 
These properties are consistent with a central star which was luminous in the early stages of its 
post-AGB evolution and therefore of high mass. Also the high central star temperature we find 
(see § §) is consistent with this. An image of the nebula is shown in Figure |l| We have obtained 
an SWS spectrum of this nebula, which is remarkable in the fact that it shows virtually no dust 
emission (see Fig. |^. This is indicative of a low F(IR)/F(II/3) ratio and dust temperature for this 
object, which confirms the evolved status of this PN. It also makes the object very suitable to 
investigate the points mentioned earlier. 

First we will discuss the observations gathered for this source in § |2|. Since the nebula is very 
large, aperture corrections are needed. These are discussed in § |3|. Next we discuss the reddening 
correction of the spectra in § ^ Some details concerning the observed spectrum, including the 
possible influence of a recombining halo, are discussed in § ||. Nebular conditions derived from 
plasma diagnostics and the possible presence of t^-effects are discussed in § y. Then we study 
the observed thermal grain emission in § [^ Finally we construct a photo-ionization model of the 
nebula in § |^, which is used to study the possible presence of dust grains inside the ionized region 
of the nebula. Final conclusions are given in § 0. 



2. Observations 

2.1. Unpublished Data 

In the analysis of NGC 6445 several spectra were used, covering a wide range of wavelengths. 
A log of all the observations that were not previously published can be found in Table |l[ The 
International Ultraviolet Explorer {lUE) spectra were taken from the INES catalog in the lUE 



final archive. No further processing was done on these data. For further reference the SWP17030 
spectrum will be called lUEl and the SWP38267 spectrum IUE2. The LWP17434 spectrum does 
not contain any detected lines and will not be discussed further. 

The optical observations were done with the Boiler &; Chivens spectrograph at the 1.52-m 
ESO telescope on La Silla. CCD #24, which is a Ford Aerospace 2048L, and grating #25 were 
used. At a grating angle of 8° this gave us a wavelength coverage from 391.4 nm to 967.8 nm, or 
0.283 nm per pixel and a resolution of about 0.7 nm. We were not affected by the second order 
spectrum with this setting. A special technique was used to obtain large-aperture exposures. For 
this the slit of the spectrograph was oriented North-South and the tracking rate of the telescope in 
the East- West direction was offset by 0.05 arcsec/s. This caused the slit to gradually drift over the 
nebula, hence causing the spectrograph to integrate photons from the entire nebula. The exposure 
time listed in Table || is the effective exposure time, i.e. the slit width divided by the differential 
tracking rate times the number of exposures. This technique allows the optical spectrum to be 
directly comparable to other large-aperture data. 

The reduction was done using standard procedures in IRAF. First the spectra were bias 
subtracted and flat-fielded with a normalized dome flat. The spectra were corrected for slit 
illumination. The images were geometrically corrected to rectify the spatial axis before sky 
subtraction. The dispersion axis is parallel to the chip's axis well within 2 pixels along the slit. 
This was good enough for our purposes as the spatial extent for the untraceable PN was much 
larger than this. The standard stars were traced when extracted. The PN spectra were extracted 
with a fixed aperture based on the emission line with the largest spatial extent. The spectra were 
corrected for airmass and sensitivity with the response curve based on the standard star LTT 7379. 
Finally the three spectra were averaged with min-max rejection. 

The infrared data were all taken with the Infrared Space Observatory (ISO), using the Short 
Wavelength Spectrometer (SWS) and the Long Wavelength Spectrometer (LWS). The data 
comprise two full grating scans (using the observing templates SWSOl with scan rate 2 and LWSOl 
with sampling interval 1) and a set of deeper exposures of selected individual lines (using the 
SWS02 template). A detailed description of the observing templates can be found in de Graauw 
et al. (1996) for SWS and Clegg et al. (1996) for LWS. The SWSOl and LWSOl data were reduced 
from the Auto Analysis Result produced with the ISO pipeline, version 6.0. The reduction was 
done using ISAP vl.5 and LIA v7.1, and consisted of the removal of cosmic ray hits and other 
artifacts and the subsequent averaging and smoothing of the data in each band. The SWS02 data 
were reduced from Standard Processed Data produced with the ISO pipeline, version 6.0. The 
SWS Interactive- Analysis System was used to remove spikes attributed to ionized-particle impacts 
and to align the continua. 

Almost all of the emission lines identified in the SWS and LWS spectra have never been 
observed before and are a valuable extension to the UV and optical spectra of this nebula. This 
is especially the case since various lines of highly ionized atoms have been detected, which are 



impossible to observe from the ground. These hnes allow a more reliable determination of the 
stellar temperature than was hitherto possible. To emphasize this point we show the SWS02 
detection of the [Ne VI] 7.65 fim line in Figure ^ which is now the highest excitation line ever 
observed in the spectrum of NGC 6445. 



2.2. Literature Data 

The optical spectrum was supplemented with data in the range from 372.6 nm to 397.0 nm 
taken from Aller et al. (1973). These data were obtained at three different positions in the nebula 
and subsequently averaged, so that they are a fair approximation for the integrated flux. The 
[Ne V] 342.6 nm flux was taken from Rowlands et al. (1993). These observations used a 48" 
aperture, which should be large enough to encompass nearly all of the [Ne V] flux. The adopted 
fluxes are listed in Table ^. 

An Ha+[N II] image of NGC 6445 was obtained by Schwarz, Corradi Sz Melnick (1992) and is 
shown in Figure |l|. From this image it is clear that the ionized region of NGC 6445 comprises not 
only the high surface brightness region that is readily observable, but also low surface brightness 
regions that extend much farther. The entire ionized region measures at least 2'x3'. 

The nebula has been observed at radio frequencies on several occasions, both using single-dish 
telescopes and interferometers. There is a significant difference between the fluxes obtained with 
the two types of telescopes. This is most likely due to the presence of extended low surface 
brightness radio emission which was included in the beam of single dish telescopes but must have 
been missed in the interferometer observations. The discrepancy between the fluxes cannot be 
attributed to optical depth effects since modeling of the nebula shows that optical depth effects 
only set in below 1 GHz. Hence it is safe to assume that the radio emission is still optically thin, 
even at the lowest frequencies observed with the VLA and Parkes. For the total radio flux emitted 
by the nebula (including low surface brightness regions) we adopt the 14.7 GHz radio flux of 
322 mJy obtained with the Parkes 64-m radio telescope (Milne & Aller 1982). 

The distance of NGC 6445 is not accurately known. Sabbadin (1986a) lists two extinction 
distances of 2.0 kpc and 2.5 kpc. Various other determinations have been published in the 
literature which usually are lower. We will adopt the statistical distance based on an empirical 
relation between the surface brightness temperature and the radius of the nebula discussed by 
Van de Steene & Zijlstra (1995). In order to obtain a distance with this method, we will need an 
angular diameter and the radio flux contained within that radius. The Parkes data are not suitable 
for this purpose, since it contains emission from low surface brightness regions and the angular 
diameter of that region is not well defined. As an alternative we will use the VLA 1.46 GHz data 
presented in Phillips & Mampaso (1988). They obtained a radio flux of 243 mJy and an angular 
extent of the high surface brightness region of 30'.' 7x40'.' 2. We will adopt the resulting geometric 
mean of 35'.' 1 as the angular diameter. It is consistent to use these data and the resulting value 



for the distance is 1.5 kpc. This result is not very sensitive to the adopted radio flux density Su 
since it scales as S~^'^^. Using the Parkes flux instead of the VLA flux would yield a result that is 
only 15% lower. 

Using the data presented above we find that the dynamical age tdyn of the nebula is 
^dyn = ©r-D/fexp = 3300 yr. Here 0r stands for the angular radius of the nebula, D for the 
distance and Uexp for the expansion velocity. 



3. Aperture Corrections 

Although all of the spectral data we present have been taken with large apertures, none of 
them encompass all of the nebula. This is especially the case for the lUE and SWS observations 
where significant parts of the high surface brightness regions have been missed. Therefore aperture 
corrections will be necessary to use these data for modeling the nebula. The apertures are 
over-plotted on an image of NGC 6445 in Figure ||. 

The lUE-SWV large aperture measures 9'.'07x21'.'65 (lUEl & 2). SWS Data in the range 
between 2.38 fiva — 12.0 ;um were taken with an aperture of 14"x20" (SWS-A), between 12.0 ;um - 
27.6 ^m with an aperture of 14"x27" (SWS-B), and between 29.5 //m - 45.2 /xm with an aperture 
of 20"x33" (SWS-C). Data in the range between 27.6 //m - 29.5 ^m were taken with a 20"x27" 
aperture. However, since there are no observed lines in this range, this aperture will be ignored in 
the subsequent discussion. To obtain aperture correction factors for the lUE and SWS data, the 
total flux was integrated inside each of the aperture boxes and compared to the total flux emitted 
by the nebula. For this we used the Ha+[N II] image of Schwarz et al. (1992). We obtained the 
following aperture correction factors which we will adopt in the remainder of the paper. For lUEl: 
15.43, for IUE2: 9.98, for SWS-A: 5.27, for SWS-B: 4.65 and for SWS-C: 3.45. 

The correction factors for the lUE data are large since this instrument measured only a small 
region in the center of the nebula, missing most of the highest surface brightness regions. This 
correction could have been overestimated due to any of the following three effects. First, errors 
of a few arcsecond in the satellite pointing could cause more flux to enter the aperture. Second, 
light from the brightest regions just outside the aperture could have been scattered into the line 
of sight by the intervening ISM. Third, the lines lUE detected come from highly ionized species 
which could be more centrally concentrated than the H I+[N II] emission. The lUE fluxes should 
therefore be considered more uncertain than any of the other data. 

The fact that the SWS apertures are centered on a surface brightness peak assures that 
the first two effects will have less influence on the correction factors for these data. Ionization 
stratification could however still be important. To study this further, we inspected the [O III] 
image taken by Schwarz et al. (1992). The correction factors for the SWS data derived from this 
image agree with the ones listed above at the 10% level. This can be explained by the fact that 
the positioning of the apertures (which include both ionized and neutral material) favors a small 



effect on the correction factors. We therefore conclude that the aperture corrections for the SWS 
data are fairly reliable. 

The apertures for the optical and ISO LWS observations are much bigger than the lUE and 
ISO SWS apertures, and cover all of the high surface brightness regions of the nebula. We will not 
apply any aperture corrections to these data. 



4. Interstellar Reddening 

The values for the interstellar extinction towards NGC 6445 quoted in the literature vary 
considerably and should therefore be considered uncertain. In view of this fact and also because 
we have spectra covering a wide range of wavelengths, a careful analysis of the data is warranted in 
order to obtain an accurate extinction curve. Such an analysis is necessarily an iterative process: 
in order to obtain an extinction curve we need estimates for the electron temperature and density, 
which in turn can only be determined accurately after an extinction correction has been applied. 
Luckily the dependence on these quantities is only moderate to weak and simple estimates are 
sufficient. In the following paragraphs we adopt values for the electron temperature and density 
which do not agree precisely with each other or with the final values which will be derived further 
on in the paper. We are however convinced that the assumed values are sufficiently accurate to 
ensure that the errors introduced by these assumptions are far less than any other source of error. 

We will start the analysis by comparing the radio and H/3 fluxes. Since we estimate that the 
optical aperture encompasses roughly 90% of the total flux emitted by the nebula, it is appropriate 
to use the Parkes 14.7 GHz flux of 322ibl8 mJy (Milne & Aller 1982) for this comparison. 
Substituting this value in Equation IV-26 of Pottasch (1984), and adopting Y = 1.23 it 0.05 and 
an average electron temperature of Te = 13.8 it 1.3 kK (obtained from a preliminary analysis 
of the [O III] lines), we were able to obtain an estimate for the dereddened (intrinsic) H/3 flux: 
log[F(iJ/3)/(Wm~^)] = —13.06 it 0.04. When we compare this with the observed H/3 flux of 
log[F(i//3)/(Wm~2)] ^ -14.12 ± 0.04, we find that c(H/3) = 1.06 ± 0.06 dex, or alternatively 
A(486.1 nm) = 2.64 mag. 

We will now determine the shape of the extinction curve for our spectra. The relative 
line strengths in the recombination spectrum of H I and He II depend very little on electron 
temperature and density and can be determined using Case B theory. This allows us to determine 
the absolute extinction at various other wavelengths using the value for A(486.1 nm). We adopt 
the intensity ratios calculated by Storey &: Hummer (1995), assuming an electron density of 
10^ cm~^ and an electron temperature of 12.5 kK and 15 kK for the H I and He II emitting regions 
respectively. These assumptions introduce an error of approximately 1% or less in the Case B line 
ratios and are therefore well justified. The results are shown in Table |2[ The Case B ratios are 
normalized to H I 486.1 nm and He II 468.6 nm for the hydrogen and helium spectra respectively. 
In order to obtain a value for the extinction at 468.6 nm, a fit was made to the observed extinctions 



for the hydrogen lines only, using a standard ISM extinction law (Cardelli, Clayton & Mathis 
1989), and adopting a ratio of the total to specific extinction Ry = A{V)/E{B — V) = 3.3. This 
allowed us to obtain an absolute flux of the He II 468.6 nm line and include the other He II lines 
in the analysis. The results are shown in Table ^ and Figure Q. 

To obtain an extinction curve, we made a least-squares fit of the law given by Cardelli 
et al. (1989) to these data, using the ratio of total to selective extinction, Ry, and the selective 
extinction, E{B — V), as free parameters. The result was Ry = 3.55 and E(B — V) = 0.61, 
which is equivalent to c(H/?) = 0.99 dex. We have composed an extinction law using the results 
of Cardelli et al. (1989) with Ry = 3.55 for wavelengths between 0.1 fira and 3.33 fj,ra, Rieke 
&; Lebofsky (1985) between 3.33 //m and 13 //m, and Mathis (1990) for wavelengths longward 
of 13 ^m. The extinction law that is obtained this way is plotted in Figure ^ and was used to 
deredden all observed fluxes. 

Altering the value of Ry has a strong effect on the UV part of the extinction curve, but 
leaves the optical and infrared part virtually unchanged. Hence the fact that we use an Ry which 
deviates from the canonical value Ry = 3.1 is important for the dereddening of the UV line fluxes 
and assures that the Case B line ratios for He II are obeyed. On the other hand, if we had used 
an extinction law with Ry = 3.1, this would only have made a difference of 10% or less for the 
optical and infrared line intensities. In other words, in cases where only optical and/or infrared 
data are available, the value of Ry cannot be determined from these data and assuming Ry = 3.1 
is appropriate. But in cases where UV data is included, Ry should be treated as a free parameter 
to assure a proper dereddening. 

The results of the dereddening can be found in Table |3| where we list the fluxes and the 
dereddened intensities relative to H/3 for all the lines used in this paper. Aperture correction 
factors have been applied to the dereddened intensities, but not to the observed fluxes (except for 
the combined lUE spectra) . The error bars quoted in Table |3| reflect purely the uncertainty in the 
measured flux and are therefore a measure for the quality of the spectrum. They do not contain 
contributions for the uncertainty in the absolute flux calibration, aperture correction or reddening 
correction. However, estimates for these uncertainties where quadratically added to make the 
photo-ionization model discussed in § |8|. 



5. Remarks on the Observed Spectrum 

The observed flux for the [O III] 88.2 /im line is more than a factor of 5 higher than the 3a 
upper limit of 2x10^^^ Wm^^ that was previously determined by Dinerstein, Lester & Werner 
(1985). This large discrepancy cannot be explained by aperture effects since both observations 
encompass essentially all of the nebula. It also seems very unlikely that it is caused by variability 
in the spectrum since there is no evidence that would support this assumption. Comparison of our 
[O III] 500.7 nm flux with the ones pubhshed in O'Deh (1963) and Aller et al. (1973) shows that 
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they all are in good agreement. Our flux is only 0.06 dex higher than the older measurements, 
making it very unlikely that the infrared [O III] spectrum has undergone large changes over the 
past 35 years. An explanation for this conundrum is therefore not apparent to us. 

The intensity /(Ha) = 363 is surprisingly large and is not expected by Case B theory. A 
possible explanation could be provided by the fact that the overall shape of the extinction curve 
can differ from the mean law for different lines of sight (Cardelli et al. 1989). HST observations 
could be used to investigate the extinction law further. 

Another explanation could be provided by the fact that the outer halo of NGC 6445 
is currently recombining (Tylenda 1986), which would cause the plasma to have an electron 
temperature well below what is expected for an ionized region and would alter the Case B line 
ratios. From our discussion in § ^ it will become clear that the central star is currently evolving 
on the cooling track in the Hertzsprung-Russell diagram, which implies that it has undergone a 
considerable drop in luminosity in recent history. This evolution could have caused the ionized 
region to shrink. Therefore it is possible that the faint surface brightness regions are actually the 
afterglow of regions that did not have sufficient time to recombine. The recombination timescale 
for hydrogen is approximately 1230T4'*^n^ yr, where T4 is the electron temperature in units of 
10^ K and 712 is the electron density in units of 10^ cm^^ (using Storey & Hummer 1995). This 
timescale could very well be longer than the evolutionary timescale of the central star, depending 
on the density in the outer halo. However, the timescale for recombination from O^"*" to O"*" is 
only 86r4'^^n^ yr (using the atomic data referenced in cloudy 90.05). This would imply that 
in the outermost regions of the halo, where the gas supposedly had the longest time to recombine, 
one would not expect [O III] 500.7 nm emission to be visible. Inspecting the [O III] image taken 
by Schwarz et al. (1992) however reveals that such emissions are present. The outermost regions 
in the North-East and South- West directions are clearly visible in both the Ha+[N II] and the 
[O III] image, indicating that they are still being photo-ionized. Some regions in the South-East 
and North- West directions however do appear much fainter in [O III], and it is plausible that 
recombination is going on there. To illustrate this point, we show in Figure ^ the ratio of the [O III] 
over the Ha+[N II] image. It can be interpreted as a measure for the level of excitation, where 
white indicates high excitation. The scale and positioning of the image is identical to Figure |l|. In 
order to increase the signal-to-noise ratio in the low surface brightness regions, the original images 
were smoothed with a gaussian with a FWHM of 2 pixels before they were divided. The resulting 
image clearly shows the complex structure of the nebula. The darkest regions in the image are 
either low excitation regions that are close to an ionization front or unilluminated regions where 
recombination is taking place. The image clearly shows that especially the North-East part of the 
halo is still highly excited. Hence the simple view where the entire halo is recombining cannot be 
correct. Due to the evolutionary status of the central star it is likely however that recombination 
is taking place to some degree. This is the case for the South-East and North- West regions, 
but probably near every ionization front as well. We cannot obtain an accurate estimate for 
the emission measure of these regions, but it seems likely it is insufficient to have a large effect 
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on the integrated hydrogen spectrum, at least within the optical aperture. On the other hand, 
recombination is expected to have a profound effect on the conditions in the Photo-Dissociation 
Region (PDR). This is because the region which is now the FDR, was in recent history part of 
the ionized region and therefore still has a 'memory' of those conditions. In particular, its energy 
content will be higher than expected from equilibrium calculations, and therefore also the electron 
temperature will be higher. This implies that cloudy may not be able to correctly predict the 
PDR emissions for this nebula. 

The fact that ionization fronts can be seen in virtually every direction in Figure |5| is a strong 
support for the assumption that the nebula is ionization bounded and that virtually no Lyman 
continuum photons can escape it. That the nebula is ionization bounded is possible because the 
central star is on the cooling track and emits far less ionizing photons than when it was still on 
the horizontal track. The fact that ionization fronts can clearly be seen in the outer halo proves 
that the halo cannot be a reflection nebula. 



6. Plasma Diagnostics and the t^-efFect 

Accurate abundance determinations are hampered by the so-called t^-problem (Peimbert 
1967). This problem is probably one of the most important in present day PN research and 
subject to much debate. The f^-parameter is a measure for the electron temperature fluctuations 
in the ionized region of the nebula. Models do not predict such fluctuations, but the limited 
amount of observational material gathered so far seems to contradict this. The i^-problem can 
be investigated by comparing the strength of infrared lines with optical lines. The emissivities of 
the far-infrared fine-structure lines have a very different dependence on electron temperature and 
density as compared with the optical lines (Simpson 1975; Watson & Storey 1980). The infrared 
line-strengths are virtually insensitive to Tg, because of the very small energy of their upper levels, 
whereas the optical lines originate from higher energy levels and have a strong dependence on Te. 
Fluctuations in Te increase the strength of emission lines and this increase will be larger for lines 
originating from higher energy levels (e.g., for optical lines). Hence the t^-effect can be measured 
by comparing electron temperatures determined from optical and infrared diagnostic lines. 

Using the dereddened line intensities, we are able to derive values for the electron temperature 
and density using various diagnostic ratios. The results are listed in Table Q. The results for the 
optical density diagnostics ([S II] and [O II]) were obtained with the package nebular in iraf 
(Shaw & Dufour 1995). All other diagnostics were obtained using our own five level model atom, 
with collision strengths taken from Lennon & Burke (1994) for [N II], [O III] and [Ne V], Butler 
& Zeippen (1994) for [Ne III], and Galavis, Mendoza & Zeippen (1995) for [S III]. The transition 
probabilities were taken from Galavis, Mendoza & Zeippen (1997) for [N II], [O III], [Ne III] and 
[Ne V], and Mendoza & Zeippen (1982) for [S III]. 

The results for the various temperature diagnostics show quite a large spread, which can 
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in part be explained by temperature stratification in the nebula. However, one also notices 
that the infrared temperature diagnostics seem to give systematically lower results than the 
optical diagnostics. This could point to a t^-effect, although the evidence is not convincing. The 
temperatures we derived indicate the following values: t^[0 III] = 0.065±0.053 and t^[S III] = 
0.168ib0.148. Hence the error margins are too large to validate the presence of a t^-effect. 

Other explanations for the difference between the optical and infrared diagnostics are also 
plausible. First of all, in order to calculate the temperature from the infrared diagnostics, the 
optical and ISO spectra need to be compared. Hence systematic effects may be introduced by 
errors in the absolute calibration of either of these spectra and/or by errors in the aperture 
corrections for the SWS data. The uncertainty introduced by the absolute calibrations has been 
estimated at 20%, and the uncertainty introduced by the aperture corrections at 15%. These 
numbers were used to calculate the error margins for the electron temperatures given above, but 
they may have been underestimated. 

Another possible explanation is inaccuracies in the atomic data, in particular the collisional 
cross sections. Questions concerning the [Ne V] data calculated by Lennon & Burke (1994) have 
been raised by Oliva, Pasquali & Reconditi (1996). However, their argument can be refuted as 
it is based on an inaccurate flux for the [Ne V] 14.3 /um line, as is shown by van Hoof et al. 
(1997b). Subsequent work (van Hoof et al. 1999) has shown that, at least for [Ne V], deviations 
up to a maximum of 30% in the data for the transitions in the ground term could be indicated by 
observations of PNe. For the data in this paper, the difference between the optical and infrared 
diagnostics can be eliminated by adjusting the collision strengths of the transitions in the ground 
term at the 30% level or less. This argument indicates that it is at least conceivable that the 
differences between the optical and infrared diagnostics are caused by inaccuracies in the collisional 
cross sections. 

The various density diagnostics are in excellent agreement, which suggests that the assumption 
of constant density is a good approximation for modeling the ionized region of this nebula. The 
density we find is fairly low for planetary nebulae, which indicates that NGC 6445 is well evolved. 



7. Thermal Grain Emission 

The fact that hardly any grain emission is detected in the SWSOl spectrum raises the question 
whether the inner regions of NGC 6445 could be free of grains, either due to grain destruction, or 
to dust-gas separation. The LWS spectrum covers a larger part of the nebula, and grain emission 
is observed in this spectrum, which indicates that at least part of the grains survived. LWS will 
have detected the emission from the grains closest to the central star, but still may have missed 
the emission from the cooler grains farther out in the nebula. Before we can study this spectrum, 
it needs to be corrected for the contribution from the ISM background. 
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7.1. Correcting the On-Source Spectrum 

We took an off-source LWS spectrum 3' north of the nebula at a position angle of roughly —50° 
with respect to the polar axis of the nebula (using the standard convention where positive angles 
indicate a counter-clockwise direction). The approximate position of the polar axis is indicated in 
Figure ||. The off-source spectrum was intended to measure the background continuum and line 
emission from the ISM. The off-source spectrum is shown in Figure g. It shows two clearly distinct 
emission components: one very cool component which only becomes detectable at wavelengths 
longward of 140 /xm, and a much warmer component which is visible at shorter wavelengths. 
The latter is observable as a smooth background in HIRES images made from IRAS 60 ^m and 
100 /im data. Measurements of this background are shown with a diamond in Figure |^. The flux 
densities were calculated from the measured surface brightness assuming a circular aperture with a 
diameter of 84" for LWS and applying the appropriate AC-DC and color corrections. The excellent 
agreement of the 100 ^m point clearly indicates that the warm component of the background 
emission is associated with cirrus. No significant gradient of the cirrus emission is seen in the 
vicinity of NGC 6445. The very cold component is not detected in the 100 ;um HIRES image since 
its flux is too low compared to the cirrus emission. The cold component is clearly extended (it is 
seen in both the on-source and off-source spectrum). Estimates for the dust temperature of both 
components are derived below. The quality of the off-source spectrum was poor at wavelengths 
shortward of 85 ^m and therefore an extrapolation was made using the flt discussed below. Hence 
the extrapolation is primarily based on the background flux measured in the 60 ^um HIRES image. 
The ISM contribution is very small compared to the on-source flux at 45 /xm (less than 2%). 
For longer wavelengths the ISM contribution becomes increasingly important: at 85 /im the ISM 
contributes 25% to the total on-source flux, at 145 /xm 50%, and at 190 ^m 90%. 

We need to correct the on-source spectrum for the contribution by the ISM background since 
this emission would lead to a serious underestimation of the nebular grain temperature. This 
correction is easily done assuming that the contribution by the ISM is the same in the on-source 
and off-source spectrum. This assumption is justifled, both for the warm and cold component of 
the ISM emission. From the HIRES images we know that the cirrus emission shows no gradient 
in the vicinity of the source. The very cold ISM component is the dominant feature at long 
wavelengths in the on-source spectrum and this feature is very nicely removed after subtraction. 
The remaining question is if grain emission from the outermost regions of the nebula itself could 
have contributed to the off-source spectrum. This argument leads to estimates of hundreds or 
even thousands of solar masses for the total nebular mass. This is clearly impossible unless we 
have seriously overestimated the distance to the nebula. This argument seems very contrived 
and we therefore dismiss it. The on-source spectrum was corrected for the contribution of the 
ISM continuum using the entire spectrum shown in Figure |6|. The corrected spectrum is shown 
in Figure 0. In the off-source spectrum two emission lines were detected, [O I] 63 /xm and [C II] 
157 /xm. The fluxes of these lines were subtracted from the on-source flux to constrain our 
photo-ionization model of the nebula. 
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7.2. Temperature Determination 

To determine the average grain temperature of the ISM background emission and the nebular 
grain emission, we will use an emission law of the form Fi, = cz/'^^i?;/ (Tdust) where c, /3 and Tjust are 
free parameters. Such a simple power law dependence of the absorption cross section is assumed 
for astronomical (i.e., amorphous) silicates beyond the 18 //m feature and for graphite beyond 
35 //m (e.g., Draine & Lee 1984). However, recent ISO observations have shown that crystalline 
silicates show a rich spectrum of emission features, both in the SWS and LWS wavelength range 
(e.g.. Waters et al. 1996), and a power law approximation is not valid for such grains. In the 
on-source spectrum we do not detect any crystalline silicate features in the LWS range, but we 
cannot rule out the presence of such features below the detection limit in the SWS observations. In 
the absence of better data, we will assume that the power law approximation is valid in the LWS 
range. This assumption is commonly made for fitting diffuse ISM emissions, often using [3 = 2. 



7.2.1. The ISM Background 

Before we discuss the grain emission from the nebula itself, we discuss the fits to the ISM 
background emission. We derived estimates for the temperature of the warm and very cold 
component by using a two-component model. Each component has the same analytic dependence 
given above and we assumed (3 = 2 for both components. To constrain the fit we used the 60 jiva 
point from the HIRES image and the LWS off-source spectrum longward of 85 ^um. The optimal 
fit is shown in Figure |y and Table p. We find that the warm cirrus component has a temperature 
of 24 K, which makes it possibly the warmest cirrus cloud known. Lagache et al. (1998) show 
various measurements in their Figure 2, the highest being 23 K in the Orion complex. The 
galactic coordinates for our measurement are ^n = 8.119, hn = -1-3.930. Comparison with Figure 1 
of Lagache et al. (1998) shows that they are in the middle of a high-temperature region. Since 
their data is smoothed over a 7° beam, small clouds with higher temperatures could have been 
averaged out. Taken at face value, our detection of a very cold component of 7 K is in agreement 
the results of Wright et al. (1991) and Reach et al. (1995), who derived the widespread presence 
of such a component based on FIRAS data taken by COBE. Lagache et al. (1998) subsequently 
dismissed these claims, upon a re-interpretation of DIRBE and FIRAS data. Since our observation 
encompasses only a small region of the sky, it is not clear how extended the 7 K component is and 
whether it could have contributed significantly to the FIRAS spectrum. Note also that our line of 
sight is in the galactic plane. Given this fact, the very cold emission component may come from 
a starless core. Ward-Thompson et al. (1994) have observed two such cores with temperatures 
of 8 K. Bernard et al. (1992) predict that for a core to reach a temperature as low as 7 K, the 
optical depth should at least be Ay = 100 mag. So for NGC 6445 to be observable, the core 
should reside behind the PN, giving a minimum distance of 1.5 kpc. Because we see the very 
cold component in both the on-source and off-source spectrum, the angular extent must be at 
least 5'. This leads to a lower limit of 2 pc for the diameter of the core. This is considerably 
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larger than the cores observed by Ward-Thompson et al. (1994). It is not clear whether cores this 
large can sustain a temperature of 7 K (see e.g. the discussion in Laureijs 1999, and references 
therein). The ISOPHOT Serendipity Survey (ISOSS) also finds very cold cores with temperatures 
of approximately 13 K and a radial extent consistent with our observations (Hotzel et al. 1999). It 
should be noted that ISOSS can only detect cores with temperatures in the range 10 K < Tjust < 
70 K (Toth et al. 1999). Hence the non-detection of cores with temperatures below 10 K in ISOSS 
is expected. Further observations are needed to clarify the nature of the very cold component. 



1.2.2. The Nebular Grains 

The LWS data for NGC 6445 show evidence for the presence of dust which is cooler than 
observed in most other planetary nebulae. We will investigate this emission first, before we try to 
determine where the emitting grains reside. The resulting best fit gave /? = 0.89 and is shown in 
Table |5| and Figure 0. Fits with /3 fixed to other values are also shown for comparison in Table g. 
The total flux emitted by the grains, F(IR) = 2.38x10"^'^ Wm~^, is nearly a factor of two lower 
than was previously reported by Lenzuni, Natta & Panagia (1989) and Zhang & Kwok (1993). 
This discrepancy can be attributed to the fact that both papers based their determination on the 
IRAS fluxes in all four bands. However, the IRAS 12 //m and 25 /xm detections of NGC 6445 are 
strongly dominated by line emission, which leads to an overestimation of the total infrared flux. 

The temperature of 57 K for the nebular dust is considerably lower than what is typically 
observed in planetary nebulae. Dust grains with such typical temperatures (roughly between 90 K 
and 150 K) cannot be present in NGC 6445, since they would be readily observable in the SWS 
spectrum. The ratio of the thermal grain emission to the H /3 flux of 27.0 is also very low. The 
presence of strong neutral lines in the spectrum indicates that the nebula is ionization bounded 
(see also § g). This argument is further strengthened by the fact that most determinations of 
the hydrogen and helium Zanstra temperatures for the central star are in good agreement (see 
Table |^. This is a good indication that the hydrogen Lyman continuum is optically thick. We 
therefore assume that no signiflcant fraction of the Lyman continuum photons escape the nebula. 
Our photo-ionization model (see § ^) predicts F{Ijy a) / F(H P) = 26.6, which implies that all 
of the grain heating could be accounted for by Lya photons. Apparently Lyman continuum 
photons contribute little to grain heating and this suggests either that the ionized region is free 
of grains, or that the ionization parameter is low, resulting in a low optical depth of the grains in 
the Lyman continuum (Bottorff et al. 1998). To investigate this further we will first construct a 
photo-ionization model of the ionized region of NGC 6445. 
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8. Photo-ionization Modeling 

We constructed a photo-ionization model of the ionized region using a modified version of 
CLOUDY 90.05 (Ferland et al. 1998) and fohowing the method described in van Hoof & Van de 
Steene (1999) with the fohowing aherations. Instead of the blackbody assumption we used model 
atmospheres calculated by Ranch (1997). Furthermore, we used a constant hydrogen density 
throughout the ionized and neutral regions of the nebula. In view of the complex structure of 
the nebula we did not constrain the model with an angular diameter, but alternatively fixed the 
density at the value determined from plasma diagnostics. We also used the electron temperatures 
determined from various plasma diagnostics as additional constraints on the model. We assumed 
no dust grains to be present in the ionized region and did not constrain the model with infrared 
continuum fluxes. We omitted the infrared fine-structure lines of [O I] from the modeling input 
since preliminary models showed a large discrepancy between the modeled and observed flux of 
these lines. The observed fluxes are approximately a factor of 3 to 5 stronger than the model 
fluxes and we feared that such a large discrepancy would upset the abundance determination. A 
possible explanation for this discrepancy could be provided by the fact that parts of the nebula 
are recombining, which would have a strong influence on the PDR physics. Cloudy is not capable 
of modeling such non-equilibrium conditions, and an underestimation of the flux of the [O I] 
fine-structure lines is expected because cloudy will underestimate the electron temperature in 
the recombining regions. See also the discussion in § H. 

An alternative explanation could be the fact that the collision strengths for the [O I] 
fine-structure lines are poorly determined because they are difficult to calculate (Bhatia & Kastner 
1995). For instance, cross sections for O'^+e collisions have recently been revised downwards by a 
factor of 5 to 50 (Bell et al. 1998), based on calculations only valid below 3000 K. Cross sections 
for O^+H and O'^+He collisions are only available for temperatures below 1000 K (respectively 
Launay & Roueff 1977 and Monteiro & Flower 1987). Since the [O I] fine-structure lines are 
mainly formed in regions with electron temperatures between 4000 K and 10 000 K, significant 
extrapolations of these (possibly uncertain) data are needed. This makes modeling of the [O I] 
fine-structure lines very difficult and may very well explain the discrepancy between the observed 
and calculated values. 

The parameters for the photo-ionization model we derived are listed in Table ^ 



8.1. The Central Star 

The temperature of the central star of NGC 6445 has been determined many times, mainly 
using the Zanstra method. An overview of all recently published determinations can be found in 
Table 0. Most of these are based on a blackbody approximation. The average of all determinations 
using the Zanstra method, or a variation thereof, is 188ib5 kK. This is in excellent agreement with 
the value that has been derived in this work, based on atmosphere models by Ranch (1997). This 
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indicates that in the case of NGC 6445 the Zanstra method works weh. This can be understood 
given the fact that the nebula is ionization bounded and optical depth effects due to dust grains 
are negligible. Hence the basic assumptions of the Zanstra method are valid in this case and the 
main uncertainties in applying the method will stem from observational error and to a much lesser 
amount from assumptions concerning the shape of the central star spectrum (i.e., the blackbody 
assumption). 

In order to check our determination of the luminosity of the central star, we also made 
predictions for the magnitude of the central star in various photometric bands. A comparison of 
the results with observed magnitudes is shown in Table §. The largest discrepancy is with the V 
magnitude measured by Jacoby & Kaler (1989). Their flux may be too low due to the possible 
presence of clouds during the observations. Gathier & Pottasch (1988) used a narrow-band 
filter which is free of emission lines. This technique is intrinsically more accurate than using a 
broad-band filter like Johnson V. Since their measurement is in excellent agreement with the 
predictions, this confirms our determination of the central star luminosity with photo-ionization 
modeling. 

When one compares the stellar temperature and luminosity with the theoretical calculations 
of Blocker (1995a, b), one finds that they indicate a core mass between 0.625 M© and 0.696 Mq, 
or alternatively a main-sequence mass between 3 M© and 4 Mq. The evolutionary age of the 
nebula predicted by the Blocker tracks is approximately 3500 yr and 900 yr for the 0.625 Mq and 
0.696 Mq tracks respectively. This is in excellent agreement with the derived dynamical age of 
3300 yr. 



8.2. Dust Grains in the Ionized Region 

Our model shows the carbon and oxygen abundances to be about equal. In view of the fact 
that the carbon abundance may be overestimated due to the uncertainties in the lUE data, it 
seems more likely that the nebula is oxygen-rich. However, we take no final standpoint on this 
issue and leave both possibilities open. It is a well-known fact that in an oxygen-rich environment 
silicate grains can be formed. These grains predominantly contain oxides of magnesium, silicon 
and iron, but also of other element such as calcium, titanium, chromium, nickel, etc. It is a lesser 
known fact that these elements can also be depleted in a carbon-rich environment. See for instance 
Beintema et al. (1996) where it is shown that sodium, magnesium, aluminum, silicon, calcium and 
iron are depleted in the carbon-rich nebula NGC 7027. It is unclear what types of grains these 
elements are depleted into. They may either be carbides (e.g., SiC has been observed in various 
objects), or they may be locked up as impurities in more abundant amorphous carbon or graphite 
grains. Which of these grain types are present in NGC 6445 is unclear, so we will simply refer to 
them as silicon-bearing grains. 

First we will investigate whether grain destruction could have taken place in the ionized region 
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of NGC 6445. When dust grains are destroyed, the constituent material is returned to the gas 
phase, making the composition essentially solar again. This has a profound effect on the spectrum 
being emitted by the plasma. Many strong lines of elements like magnesium, aluminum, calcium, 
vanadium, chromium, etc. should be easily detectable in UV or optical spectra (Kingdon, Ferland 
(fe Feibelman 1995; Kingdon & Ferland 1997). The fact that none of these lines are detected in 
the spectrum of NGC 6445 already indicates that no substantial amount of grain destruction (at 
least for the silicon-bearing grains) can have taken place in the ionized region. This is confirmed 
by our photo-ionization modeling which clearly indicates that elements like magnesium, silicon 
and calcium are depleted in the ionized region. The derived abundance for magnesium is marked 
uncertain because it is based on only one emission line belonging to a subordinate ionization stage, 
which implies that a substantial correction is needed for unobserved ionization stages. By itself the 
magnesium under-abundance is not convincing evidence for depletion. It is however supported by 
the under-abundance for silicon which is better determined. It is also based on only one emission 
line, [Si II] 34.8 /xm, but since Si"*" is the dominant ionization stage, only a small correction for 
unobserved ionization stages is needed. Hence the dominant source of error for the abundance will 
be the uncertainty in the measured line flux and to a lesser extent the aperture correction. It is 
very unlikely that the strength of the silicon line could have been underestimated by a factor of 
five and it is even less likely that the aperture correction could be off by such a large factor. The 
calcium under-abundance is based on a 3a upper limit for the [Ca II] 729.1 nm line. Since Ca^ is 
a subordinate ionization stage, the largest source of error will be the correction for unobserved 
ionization stages. It seems highly unlikely however that this correction could be off by a factor of 
twenty. 

Recently Liu (1998) published an abundance analysis of the planetary nebula NGC 6153, 
which yielded abundances derived from recombination lines that were roughly a factor of 10 higher 
than those derived from collisionally excited forbidden lines. Taken at face value, this could cast 
doubt on the analysis presented above. Similar discrepancies have been observed in other nebulae 
as well, but usually they are much less pronounced. The strength of this effect for NGC 6445 has 
not been established. Whether such a discrepancy exists for magnesium, silicon or calcium has not 
been established either since the recombination lines of these lines are too weak to be observable. 
The physics causing the discrepancy between collisional and recombination abundances is not 
clear and therefore both values should be viewed with caution until such an explanation is at 
hand. Liu (1998) claims that there is no evidence that the recombination abundances have been 
grossly overestimated. However, the derived O^^ abundance of 9.65 for NGC 6153 is much larger 
than the solar oxygen abundance. Significant enhancement of oxygen in post-AGB stars is not 
expected by stellar evolution theory. We argue that this fact does cast doubt on the recombination 
abundances. On the other hand, the O^"*" abundance derived by Liu (1998) from forbidden lines 
(using essentially the same method we use in this paper) is in good agreement with the expected 
solar value. Exactly the same argument holds for the Ne^^ abundances of NGC 6153. In view of 
these facts it seems fair to assume that the collisional abundances derived in this work are valid. 
Therefore we claim that there is evidence for depletion in NGC 6445, which in turn proves that 
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the dust grains that contain these elements cannot have been destroyed in situ, at least not in 
significant amounts. 

An alternative assumption could be that the ionized region is free of dust grains due to the 
fact that the grains were separated from the gas. If dust-gas separation occurred, it must have 
been during the post-AGB phase. It could not have occurred during an earlier phase, since then 
dust formation would still be on-going. On the other hand, separation during the early post-AGB 
phase could be expected as it is radiation pressure on dust grains that is thought to drive the 
outflow of gas and dust at that stage. However, this phase is very short-lived and it is not likely 
that a detectable amount of dust-gas separation could occur in that period. Separation during 
later phases seems impossible since our photo-ionization model indicates that once the nebular 
material is ionized, the drift velocity with respect to the gas is lower than 1 kms~^. This velocity 
is far too low to produce any significant dust-gas separation. 

We have constructed a model for the grain emission of NGC 6445 based on the assumption 
that it contains grains at a constant dust-to-gas ratio throughout the nebula. Heating of the 
grains by Lya photons must be included in this model. Lya photons are scattered very efficiently 
by neutral hydrogen due to Rayleigh scattering. This process tends to lock up the photons inside 
the neutral envelope of the nebula until they are eventually absorbed by grain particles. However, 
Lya photons can escape provided they can pick up sufficient redshift by repeatedly being bounced 
off the expanding neutral envelope. We have included a simplified analytic treatment of this 
process in the model. This treatment assumes that there is a sharp boundary between a fully 
ionized and a fully neutral region. Hence there is no Rayleigh scattering inside the ionized region, 
only outside. It is known that this is not a good approximation and it will overestimate the escape 
probability to some extent. We also included sources of grain heating other than Lya in our 
treatment by using the Cloudy model. These calculations indicate that roughly 25% of the total 
grain heating is due to sources other than Lya and hence that roughly 25% of the Lya photons 
must escape in order to account for the observed total grain emission. This shows that Lya escape 
should be included in the model. Given that there are considerable uncertainties in the shape of 
the absorption law for the grains and in the geometry of the nebula, we felt that a more detailed 
treatment of Lya scattering (e.g., using a Monte Carlo code) was not warranted. The shape of the 
absorption law is important because the cross section at 121.6 nm (combined with the dust-to-gas 
ratio) determines the escape probability, and because a large infrared emissivity (when compared 
to the UV absorption) will tend to make the grains cooler. Also deviations from /3 = 2 in the LWS 
range will have an important effect on the modeling. None of these parameters are accurately 
known for the grains in NGC 6445. 

We produced 3 models using single-sized grains of astronomical silicate (Martin & Rouleau 
1991) with diameters 0.01 /im, 0.1 ^m and 1 fira respectively. We treated the dust-to-gas ratio as 
a free parameter, but assumed it to be constant as a function of radius. The results are shown in 
Figure ^. The dust-to-gas ratios for the best fitting models are: 3.0x10"^ (0.01 ^m), 4.1 xlO^^ 
(0.1 ^m) and 1.7x10"^ (1 /um). The models using the 0.01 /im and 1 /im grains are clearly 
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unrealistic, but the 0.1 /im grains give a reasonable fit. The sequence of models suggests that the 
grains in NGC 6445 are on average large. 

We mentioned earlier that our model will overestimate the escape probability of Lya photons 
to some extent. Since the total amount of heating of the grains is constrained by the observations, 
the model fitted to the observations will tend to overestimate the dust-to-gas ratio somewhat 
to counteract this effect. If we had assumed the silicon-bearing grains to be crystalline silicates 
instead of astronomical silicate and if we assume the UV absorption characteristics of both silicate 
types to be the same, the crystalline silicates would have been cooler due to the presence of 
strong emission features in the mid- and far-infrared. Under those conditions our model sequence 
would have led to a somewhat lower estimate for the average grain size and dust-to-gas ratio. 
Using graphite grains would have led to a model and a dust-to-gas ratio that is roughly similar 
to what is quoted above. Given the uncertainties in our model and the absorption and emission 
characteristics of the grains in NGC 6445 we cannot derive precise results. However, it is very 
likely that the grains in NGC 6445 are large and that the dust-to-gas ratio is several times 10~^. 
One can also predict the dust-to-gas ratio assuming that the silicon abundance before grain 
formation was solar and the grains are composed of pure MgSiFe04. This yields 3.2x10"^, in 
excellent agreement with our results. We find there is no need to invoke dust-gas separation in 
this nebula in order to obtain a good fit to the spectrum. 

Therefore, the most likely conclusion is that grains are residing inside the ionized region 
of NGC 6445 and that the low temperature is caused by the low luminosity of the central star 
combined with the large extent of the nebula. The low flux of the grain emission, when compared 
to H/3, is caused by the low optical depth of the grains in the Lyman continuum, which makes 
Lya photons the dominant heating source for the grains. In short, all of these characteristics are 
a consequence of the fact that the nebula is well evolved. Our conclusion implies that the bulk of 
the silicon-bearing grains in this nebula could survive exposure to hard UV photons for at least 
several thousands of years. This contradicts the results of Natta & Panagia (1981) and Pottasch 
(1987), who find a tight correlation between the dust-to-gas ratio and the nebular radius based 
on a small sample of PNe. The dust-to-gas ratio in NGC 6445 is much higher than would be 
expected from this correlation. The results given by Natta & Panagia (1981) predict a dust-to-gas 
ratio of approximately 2x10^^ for NGC 6445, i.e. roughly 2 orders of magnitude lower than what 
we find. We are confident that our models could not have overestimated the dust content of the 
nebula by such a large factor. We cannot rule out that some destruction of the smallest grains, 
constituting only a small mass fraction, could have occurred. Small grains would be the first to 
be destroyed since they absorb UV radiation more effectively and are therefore hotter than larger 
grains. Such an occurrence would cause the grains in NGC 6445 to be on average larger. It is not 
clear whether this is sufficient to explain the large average grain size discussed above, or if it is 
simply the natural result of grain growth (and coagulation) during the prior AGB phase. Natta & 
Panagia (1981) found a tight decreasing relation between the average grain size and the nebular 
radius based on the same sample mentioned earlier. This relation would predict an average grains 
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size of 0.007 fim for NGC 6445, which is significantly smaher than what we find. Hence this 
prediction also seems implausible. 



9. Conclusions 

1. We presented ISO-SWS and LWS spectra of the planetary nebula NGC 6445 which indicate 
the presence of very cool dust with a low flux compared to H/3. The fact that Lya photons 
alone have sufficient energy to account for all of the grain heating, suggests that either the 
optical depth of the grains in the Lyman continuum is very low, or that the ionized region is 
free of grains. 

2. A photo-ionization model of the ionized region of NGC 6445 was constructed. Based on this 
model we showed that it is very unlikely that either grain destruction or dust-gas separation 
has taken place in the ionized region. 

3. Hence, the most likely conclusion is that grains reside inside the ionized region of NGC 6445 
and that the low temperature and flux of the grain emission is caused by the low luminosity 
of the central star and the low optical depth of the grains. This implies that the bulk of the 
silicon-bearing grains in this nebula could survive exposure to hard UV photons for several 
thousands of years, contradicting previously published results. 

4. A comparison between optical and infrared diagnostic line ratios gave a marginal indication 
for the presence of a t^-effect. However, the evidence is not convincing and alternative 
explanations could be provided by uncertainties in the absolute flux calibration of the 
spectra, the aperture corrections that have been applied or the coUisional cross sections. 

5. The photo-ionization model allowed an accurate determination of the central star 
temperature using model atmospheres by Ranch (1997). The resulting value of 184 kK is in 
good agreement with the average of all published Zanstra temperatures. 

6. The off-source spectrum taken with LWS clearly shows the presence of a warm cirrus 
component with a temperature of 24 K as well as a very cold component with a temperature 
of 7 K. Since our observation encompasses only a small region of the sky, it is not clear 
how extended the 7 K component is and whether it contributed significantly to the FIRAS 
spectrum taken by COBE, as reported by Wright et al. (1991) and Reach et al. (1995), but 
dismissed by Lagache et al. (1998). Because our line of sight is in the galactic plane, the 
very cold component could be a starless core. 
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Table 1. Log of the Observations. 









Exposure 


Aperture^ 


PAi= 


Central Position 


Spectral Range 


Telescope 


Instrument 


Date 


(s) 


(arcsec^) 


(degree) 


(2000.0) 


(/xm) 


lUE 


SWP17030 


820524 


10800 


9.07x21.65 


156.9 


17 49 15.0 


-20 00 34 


0.115-0.198 




SWP38267 


900228 


3600 


9.07x21.65 


165.5 


17 49 14.9 


-20 00 26 


0.115-0.198 




LWP 17434 


900227 


3600 


9.91x22.51 


165.5 


17 49 14.9 


-20 00 26 


0.185 - 0.336 


ESO 1.52m 


B&Ch 


950516 


120^^ 


60x140 


0.0 


17 49 14.8 


-20 00 40 


0.391 - 0.968 


ISO 


SWSOl 


970317 


1912<* 


20x14 -33x20^= 


91.1 


17 49 14.4 


-20 00 24 


2.35 - 45 




SWS02 


970317 


1828^1 


20x14-33x20° 


91.1 


17 49 14.4 


-20 00 24 


f 




LWSOl 


970317 


640^ 


84" circ. 




17 49 14.4 


-20 00 24 


43.2 - 189 




LWSOl OFF 


970317 


640^ 


84" circ. 




17 49 15.3 


-19 57 34 


43.2 - 189 



^The dimensions are EWxNS assuming PA=0°. 

''The aperture rotates counter-clockwise for increasing position angle. 
■^This is the effective exposure time (see § El). 
''This is the Target Dedicated Time. 
''The aperture size depends on the wavelength (see § m). 
Selected emission lines were observed (see Table Bh. 



Table 2. Empirically Determined Values for the Extinction. 





A 


Obs. Flux 


CaseB 


Pred. Flux 


1/A 


A{\) 


Spectrum 


(nm) 


(10-1^ Wm-2) 


Ratio 


(10-1^ Wm-2) 


(/.m-i) 


(mag) 


HI 


954.6 


0.146 


0.0357 


0.311 


1.048 


0.82 


HI 


922.9 


0.0922 


0.0249 


0.216 


1.084 


0.93 


HI 


656.3 


5.19 


2.8158 


24.5 


1.524 


1.68 


Hell 


541.2 


0.0508 


0.0787 


0.327 


1.848 


2.02 


HI 


486.1 


0.762 


1.0000 


8.70 


2.057 


2.64 


Hell 


468.6 


0.368 


1.0000 


4.16 


2.134 


2.63 


HI 


434.0 


0.334 


0.4713 


4.10 


2.304 


2.72 


HI 


410.2 


0.115 


0.2612 


2.27 


2.438 


3.24 


HI 


397.0 


0.0762 


0.1605 


1.40 


2.519 


3.16 


Hell 


164.0 


0.399 


6.9882 


29.0 


6.098 


4.65 
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Table 3. The observed line fluxes and dereddened intensities for NGC 6445. Where 
appropriate 3 a upper limits are given. 





A 


Flux 


a 


Intensity 




Spectrum 


(nm/^m) 


(10-14 Wm-2) 


(percent) 


/(H/3) = 100 


Comment 


CIV 


154.9 


0.015 


17.1 




lUEl 






0.026 


32.1 




IUE2 






0.247:'^ 


16.2 


283.:^^ 


combined 


Hell 


164.0 


0.024 


13.7 




lUEl 






0.040 


17.5 




IUE2 






0.399:^ 


13.2 


416.:^ 


combined 


cm] 


190.9 


0.037 


5.6 




lUEl 






0.053 


8.7 




IUE2 






0.556:'^ 


5.7 


712.:^^ 


combined 


[NeV] 


342.6 


0.050 


20.0 


13.9 




[Oil] 


372.6 


0.838 




209. 




[Oil] 


372.9 


0.769 




192. 




[Ne III] 


386.9 


0.983 




231. 




[Ne III] 


396.8 


0.312 




69.8 




HI 


397.1 


0.076 


20.3 


17.0 




[SII] 


407.2 


0.055 


51.2 


11.6 




HI 


410.2 


0.115 


19.4 


23.9 




HI 


434.0 


0.334 


10.4 


60.0 




[0 III] 


436.3 


0.123 


24.6 


21.8 




He I 


447.1 


0.058 


38.7 


9.6 




Hell 


468.7 


0.368 


6.2 


53.3 




[Ar IV] 


471.3 


0.044 


36.2 


6.2 




[Ar IV] 


474.1 


0.030 


59.7 


4.2 




HI 


486.1 


0.762 


5.0 


100. 




[0 III] 


495.9 


3.84 


4.6 


479. 




[0 III] 


500.7 


12.2 


6.5 


1490. 




[NI] 


520.0 


0.145 


18.0 


16.0 




Hell 


541.3 


0.051 


25.3 


5.1 




[Nil] 


575.5 


0.112 


12.8 


10.0 




He I 


587.6 


0.178 


8.2 


15.3 




[0 1] 


630.0 


0.503 


3.6 


37.9 




[S III] 


631.2 


0.037 


21.2 


2.8 




[0 1] 


636.4 


0.153 


4.6 


11.3 




[Nil] 


654.8 


3.19 


8.3 


224. 




HI 


656.3 


5.19 


6.3 


363. 
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Table 3~Continued 








A 


Flux 


a 


Intensity 




Spectrum 


(nm/^m) 


(10-14 Wm-2) 


(percent) 


/(H/?) = 100 


Comment 


[Nil] 


658.4 


9.51 


4.6 


661. 




He I 


668.0 


0.058 


14.1 


3.9 




[SII] 


671.6 


0.803 


5.0 


53.6 




[SII] 


673.1 


0.755 


5.2 


50.3 




[ArV] 


700.6 


< 0.052 




< 3.2 




He I 


706.5 


0.056 


7.7 


3.4 




[Ar III] 


713.6 


0.749 


3.1 


44.2 




[Ca II] 


729.1 


< 0.017 




< 1.0 




[Oil] 


731.9 


0.221 


9.2 


12.4 




[Oil] 


733.0 


0.185 


10.2 


10.3 




[Ar IV] 


733.2 


< 0.019 




< 1.0 




[Ar III] 


775.1 


0.171 


12.8 


8.5 




[Ca II] 


849.8 


< 0.033 




< 1.3 




[Ca II] 


854.2 


< 0.035 




< 1.4 




[SIII] 


906.9 


1.11 


3.7 


39.9 




HI 


923.2 


0.092 


21.1 


3.2 




[SIII] 


953.2 


2.73 


7.7 


91.1 




HI 


954.9 


0.146 


41.6 


4.9 




HI 


4.1 


0.094 


3.5 


7.2^ 


SWS02 


[Mg IV] 


4.5 


0.074 


8.6 


5.5^ 


SWS02 


[Ar VI] 


4.5 


< 0.051 




< 3.8^ 


SWS02 


[MgV] 


5.6 


< 0.142 




< 10.4^ 


SWS02 


[Ne VI] 


7.6 


0.052 


10.3 


3.8^ 


SWS02 


[Ar III] 


9.0 


0.339 


14.9 


27.4^ 


SWSOl 


[SIV] 


10.5 


0.954 


9.7 


77.0^ 


SWSOl 


[Ne II] 


12.8 


0.377 


1.6 


24.7^ 


SWS02 


[MgV] 


13.5 


< 0.032 




< 2.P 


SWS02 


[NeV] 


14.3 


1.64 


1.8 


105.*^ 


SWS02 


[Ne III] 


15.6 


4.05 


1.0 


259.^ 


SWS02 


[SIII] 


18.7 


0.796 


11.9 


51.5^ 


SWSOl 


[Ar III] 


21.8 


< 0.276 




< 17.7^ 


SWSOl 


[NeV] 


24.3 


1.98 


1.0 


126.'" 


SWS02 


[OIV] 


25.9 


10.8 


1.2 


686.^ 


SWS02 


[S III] 


33.5 


1.44 


14.3 


67.0^ 


SWSOl 


[Si III 


34.8 


0.488 


26.4 


22.6^ 


SWSOl 
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Table 3 — Continued 





A 




Flux 


a 


Intensity 




Spectrum 


(nm/;[ini) 


(10- 


-14 Wm-2) 


(percent) 


/(H/3) 


= 100 


Comment 


[Ne III] 


36.0 




0.360 


5.4 




16.7^ 


SWS02 


[OIII] 


51.7 




19.0 


1.8 




255. 




[N III] 


57.3 




6.30 


1.6 




84.4 




[01] 


63.2^ 




4.71 


1.4 




63.0 










0.124 


23.4 




1.7 


off source 


[OIII] 


88.2 




10.9 


1.4 




146. 




[Nil] 


121.8 




0.402 


2.9 




5.4 




[01] 


145.6^ 




0.235 


2.3 




3.1 




[CII] 


157.7 




0.894 


1.7 




12.0^= 










0.241 


14.0 




3.2^ 


off source 



^An aperture correction factor lias been applied. 

This line was not used to constrain the photo-ionization model. 
'^The off-source flux was subtracted from the on-source flux to constrain the modeling. 



Table 4. Plasma diagnostics for NGC 6445, using various line ratios. 







% 


ne 


Diagnostic 


Ratio 


(kK) 


(cm-3) 


[N II] (654.8+658.4)/575.5 


88±12 


10.4±0.6 




[S III] (906.9+953.2)/631.2 


48±10 


11.5±1.1 




[0 III] (495.9+500.7)/436.3 


90±23 


13.5±1.2 




[S III] 953.2/18.7 


1.8±0.5 


8.9±1.6 




[0 III] 500.7/51.9 


5.8±1.2 


11.5±0.9 




[Ne III] 386.9/15.6 


0.89±0.26 


11.0±0.9 




[Ne V] 342.6/14.3 


0.13±0.04 


12.0±0.8 




[S II] 671.6/673.1 


1.07±0.08 




470±170 


[0 II] 372.6/372.9 


1.09±0.11 




500±150 


[S III] 18.7/33.5 


0.77±0.14 




310±200 


[0 III] 51.9/88.2 


1.75±0.18 




550±90 
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Table 5. Parameters for the fits to the far- infrared 

emission (45 nm - 190 ^m) of the nebular grains and ISM 

grains respectively. 







Nebular j 


grains 






^dust 




FIR 




p 


(K) 


(10- 


-14 ^j^-2) 


XVN 


0.00 


84.6 




301 


2.21 


0.50 


66.5 




259 


1.59 


0.89 


57.3 




238 


1.45 


1.00 


55.2 




233 


1.46 


1.50 


47.3 




215 


1.74 


2.00 


41.5 
T 

-•■warm 


ISM backi 


201 
ground 

T'cold 


2.32 



P (K) (K) xVN 

2.00 24.4 6.7 0.92 
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Table 6. Parameters for the CLOUDY model of NGC 6445. All 
symbols have their usual meaning. 

Parameter Value rel. to solar^ 

Teff (kK) 184. 

L, (Lq) 1035. 

Tin (mpc) 79. 

rstr (mpc) 183. 

nu (cm-3) 525^^ 

Te (kK) 11.49^ 

ric (cm^'^) 571^^ 

e(He)'i 11.18 +0.19 

e(C) 8.87: +0.32: 

e(N) 8.39 +0.42 

e(0) 8.87 +0.00 

e(Ne) 8.31 +0.23 

e(Mg) 7.23: -0.35: 

e(Si) 6.86 -0.69 

e(S) 6.90 -0.43 

e(Ar) 6.58 +0.06 

e(Ca) < 5.09 < -1.27 

D (kpc) 1.5^''' 

X^ 2.9 

^Grevesse, Noels & Sauval (1996). 

This parameter was not varied. 
*^This is a weighted average over the entire nebula. 

e(^) stands for the logarithmic abundance of A; e(H) = 12. 
''Van de Steene & Zijlstra (1995). 



32 



Table 7. Various determinations of the stellar temperature. 



rz(Hi) 

(kK) 



Tz(He II) 

(kK) 



Other Methods 
(kK) 



Ref 



250. 
184. 



175. 
162. 

224. 

183. 



177. 

182. 

178. 
174. 



182. 

184. 
255. 
186. 



175. 
166. 
186. 

215. 



180. 

214. 
182. 
159. 



174. (ZM) 



145. (EB) 
180. (He) 



172. (He) 

186. (ZA) 
184. (ZC) 

194. (EB) 
182. (ZC) 

166. (EB) 



1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 
17 
17 
18 



References — 1. Harman &: Seaton 1966, 2. Pilyugin, Sakhibullin & Khromov 1978, 3. Koppen 
& Tarafdar 1978, 4. Pottasch 1981, 5. Kaler 1983 (using a covering factor ^ = 1.00.), 6. Preite- 
Martinez & Pottasch 1983, 7. Che & Koppen 1983, 8. Reay et al. 1984, 9. Tylenda 1986, 10. 
Sabbadin 1986b, 11. Gurzadyan 1988, 12. Gathier & Pottasch 1988, 13. Tylenda 1989, 14. 
Kaler & Jacoby 1989, 15. Jacoby & Kaler 1989, 16. Preite-Martinez et al. 1989, 17. Mallik 
1991, 18. Mal'kov, Golovatyj k Rokach 1995. 

Methods — EB: Energy Balance or Stoy method. He: Empirical relation between 7(468.6 nni) 
and T*, ZA: Average of hydrogen and helium Zanstra temperature, ZC: Zanstra temperature using 
cross-over technique, ZM: Helium Zanstra temperature method using atmosphere models. 
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Table 8. Predicted and observed central star magnitudes. 





U 


B 


479.3^^ 


V 


R 


I 


model 


14.82 


16.30 


16.43 


16.66 


16.85 


17.10 


.4(A) (mag) 


3.24 


2.74 


2.52 


2.14 


1.68 


1.14 


predicted 


18.06 


19.04 


18.95 


18.80 


18.53 


18.24 


observed 






18.97^ 


18.9^ 
19.04^ 







''ESO 1402 filter centered on 479.3 nm 
^'Gathier & Pottasch 1988 
^Reay et al. 1984 

'iJacoby & Kaler 1989 



^y.'*' 



1 arcmiri 



'.^':'}^m^; 










^. POLAR AXIS 



N» 



-.J 



Fig. 1. — The various apertures over-plotted on an Ha+[N II] image of NGC 6445 (Schwarz et 
al. 1992). The intensity has been scaled logarithmically to emphasize the low surface brightness 
regions. 
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Fig. 2.— The /^O-SWS spectrum of the planetary nebula NGC 6445. 
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Fig. 3.— The detection of [Ne VI] in NGC 6445. 
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Fig. 4.— The extinction law for NGC 6445. 










Fig. 5.— The ratio of the [O III] over Ha+[N II] image of NGC 6445 (Schwarz et al. 1992). 
Dark regions indicate fow excitation and are interpreted as regions near an ionization front or 
unihuminated regions that are currently recombining. 
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Fig. 6. — The off-source spectrum near NGC 6445. 



50 







d 


^ 


1 1 1 1 1 1 1 1 1 1 








o 


s 


observed spectrum 


- 






- 


o 


fitted spectrum, (5 = 0.89 

Sil. 0.01 micron 


- 


T 

a 


100 






^ Sil. 0.1 micron 

K Sil. 1 micron 

O 


- 


'a 














^ 




1 \ 








- 


T 
o 


50 


1 ^ "s ' 

- , / \ \ 








- 


^ 

fc 




- 1 / ,.-" \«j»j.y 

1 / /' 
" 1 / ■' 


t 


w.. 




- 







■;/^ .-'-^vj 




: 




1 






1 , . , , 1 . , , 


— -A^ 



100 

X(\xm) 



150 



Fig. 7. — The far-infrared spectrum of NGC 6445. 



